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ABSTRACT 

The mass and radius of the neutron star (NS) in low-mass X-ray binaries can be obtained by fitting 
the X-ray spectrum of the NS in quiescence, and the mass and radius constrains the properties of dense 
matter in NS cores. A critical ingredient for spectral fits is the composition of the NS atmosphere: 
hydrogen atmospheres are assumed in most prior work, but helium atmospheres are possible if the 
donor star is a helium white dwarf. Here we perform spectral fits to XMM, Chandra, and ROSAT 
data of a quiescent NS in the globular cluster M13. This NS has the smallest inferred radius from 
previous spectral fitting. Assuming an atmosphere composed of hydrogen, we find a significantly 
larger radius, more consistent with those from other quiescent NSs. With a helium atmosphere (an 
equally acceptable fit), we find even larger values for the radius. 

Subject headings: binaries : X-rays — dense matter — stars: neutron — globular clusters: individual 
(NGC 6205) 



1. INTRODUCTION 

Understanding the interiors of neutron stars (NSs) by 
measuring their mass es and radii is a key goal o f high- 
energy astrophysics (jLattimer fc Prakashl 12007ft . NS 
masses can be effectivel y measured for radio pulsars (e.g. 
Antoniadis et al.1 120121) and some X-ray binaries (e.g. 
Rawls et alll2011j ) by radio timing and/or optical/IR ra- 
dial velocities. However, NS radii are much more com- 
plicated to measure. The radius measured at infinity 
is affected by the NS's gravitational redshift, as Roo = 
R(l + z) = R/y/1 - 2GM/{Rc 2 ) (where R and M refer 
to the NS values), and in more subtle ways by the surface 
gravity, magnetic fields, emission anisotropics, and com- 
position. Two key methods of constr aining NS radii hav e 
involved fitting X-ray burst spectra (|Damen et al.lfT990h . 
and fitting X-ray spectra of quiesc ent low-mass X-ray b i- 
naries (qLMXBs) containing NSs (jRutledge et alJll999h . 
Recently, substantial work has been done using RXTE 
burst spectral measurements to constrain the mass and 
radius of NSs where the distance can be inferred, though 
the interpretation of these resu lts has differed due to 
different choices of assumptions ( Boutlou kos et ai]|2010t 
Ozel et al.1 120101: iSteiner et all 120101: iSuleimanov et al.1 
201lHZamfi7 et al .1 120121 : IGallowav fc Lampell2012h . 

The qLMXBs typically show thermal blackbody-likc 
radiation and/or a harder nonthermal component, of- 
ten fit with a powe r-law of photon index 1 to 2 
(|Canrpana et al.lll998|) . The nonthermal component is 
of uncertain origin, possibly due to continued accretion 
or a pulsar wind. The thermal component is better un- 
derstood, as emission from the NS surfa ce, powered by 
some combination of continued accretion (IZampier i et al.l 
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119951 ). he at release from the cru s t deposited by the las t 
outburst (jRutledge et al.ll2002bt iDegenaar et all 120111 ). 
and heat deposited in the core (created by nuclear fusion 
in the d eep crust during o utbursts) that is now leak- 
ing out (|Brown et al.lll9 98). Fitting the thermal X-ray 
spectra of qLMXBs is conceptually simpler than fitting 
X-ray bursts, as accretion has either stopped or is at a 
very low level, allowing the stratification of the lightest 
element to the top o f the a tmosphere within 30 seconds 
([Alcock fe Illarionovl I1980D . Calculations of nonmag- 
netic hydrogen atmosphere NS spectra generally reach 
very go od agreement with each other and with obser- 
vations (iRaiagopal fe Romanil Il996t iZavlin et al.l I199& 
iHeinke et al.l 120061 : lHaakonsen et al.l 120121 ). suggesting 
that X-ray spectra of qLMXBs can provide useful NS 
mass/radius constraints. 

However, qLMXB radius constraints depend strongly 
on the distance, and the distances to most qLMXBs arc 
poorly known. One solution is to study qLMXBs in glob- 
ular clusters, where t he distances are typically known to 
~ 5 - 10% precision (jBrown et all 1 199a IRutledge et all 
l2002al) . A number of qLMXBs have b een studied in 
globular clusters (see iGuillot et all 12009!. for a review), 
but only a few have sufficiently high-quality spectra 
to provide interesting constraints on the equation of 
state: these include X5 and X7 in 47 Tu c (IHeinke et aF 
20034 12006IL and the q LMXB s i n u Cen (IRutledge et af 
2002at iGendre et all I2003bt IWebb fc Barret| I2007fr 



NGC 6397 (IGrindlav et all 120011: IGuillot etHf 1201 lal) 



M28 dBecker et al.l 1 20031: IServillat et~aTl [21)121) . M13 
(IGendre et alll2003al IWebb fc Barretl 120071) . and NGC 
6553 (jGuillot et alll2011bl )~ Interestingly, two of these 
constraints (X7 in 47 Tuc, vs. M13) disagree with each 
other at the 99% confidence level, motivating efforts to 
understand the discrepancy. U ncertainties in the pileup 
correction used for 47 Tuc X7 (|Davisll20"o"lj) might drive 
the discrepancy, suggesting a deep observation of 47 Tuc 
with a smaller frame time to eliminate pileup. Alter- 
natively, the discrepancy may be due to differing atmo- 
spheres. 
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A large fraction of bright persistent or transient 
LMXBs in globular clusters have orbital periods less than 
one hour (so-called "ultracompact" systems, 5 of the 
10 systems with known periods), requiring degenerate 
white dwarf companions; such systems can be easily cre- 
ated through dynamical interactions in g lobular clusters 
()Deutsch et al.ll2000t Hvanova et al.ll2005r i. Ultracompact 
systems are likely to accrete material without hydrogen, 
since the donor star is (for the most typical evolutionary 
tracks) devoid of hydrogen, though spalla tion of infalling 
mate rial to create hydrogen is possible (jBildsten et al.l 
1 1992D . Several characteristics of X-ray bursts differ be- 
tween ultracompact vs. long-period LMXBs: the ratio of 
integrated persistent to burst fluence (much larger for ul- 
tracompact systems, indicating a lack of H) , the existence 
of intermediate-long bursts at low mass transfer rates 
in ultracompacts (ignition of large He columns without 
H bursts), higher Eddington limits from ultracompacts 
than long-period sources (where distances are indepen- 
dently estimated from their globular cluster locations), 
and consistently short bursts at moderate mass transfer 
rates in ultracompacts ( when mixed H/He ignition is ex- 
pected if H is present ) (|Cumming| 120031 : lin't Zand et al.l 
12001 iGallowav et al.l 120081 120101 ) . Thus it seems likely 
that many globular cluster qLMXBs may have atmo- 
spheres composed of He, or C, depending on the nature 
of the donor star. 

The possibility of differing atmospheric composition 
for globular cluster NS qLMXBs may explain the dis- 
crepancy between globular cluster qLMXB mass/radius 
constraints, since fits with He atm osphere models giv e 
larger radii than H atmosphere fits (|Ho fc Heinkdl2009l ). 
We used the qLMXB in M28 as a first example of 
the diffe rences in mass /radiu s constraints from the two 
models (jServillat et al.l 120121 ). Here we consider the 
qLMXB with the smallest known radius constraint, the 
M13 qLMXB. M13 has been studied by ROSAT's PSPC 
and HRI cameras, which detected an X-ray source (la- 
beled Ga) in the core of the cluster dFox et al.l 119961: 
iVerbuntl 120011 ). the target of this study. iGendre et al.1 
(|2003af ) used XMM to identify another source, 15" to 
the NW of Ga, which may contribute to the ROSAT 
PSPC and XMM spectra of Ga, and showed that Ga's 
spe ctrum was consistent with a hydrogen-atmosphere 
NS. lWebb fe Barretl (|2007l ) then calculated detailed con- 
straints on its mass and radius, along wi th XMM stud- 
ies of two other NSs in globular clusters. iServillat et al.l 
(|2011[ ) used ground-based optical (Faulkes Telescope 
North), Chandra X-ray, and Hubble Space Telescope 
data to identify the source NW of Ga (their star 4, or X6; 
we use the latter name) as a cataclysmic variable (CV) 
experiencing a dwarf nova outburst. Some additional 
results from the archival Chandra observations o f M13 
(PI: Lewin) ha ve been published (jPoolev fe Hut! 120061 : 
iHui et all 120091 ) . but spectral analysis of Chandra data 
on the NS qLMXB has not yet been published. 

2. OBSERVATIONS AND DATA REDUCTION 

Our principal dataset is a pair of XMM- Newton obser- 
vations obtained in Jan. 2002 using the EPIC cameras 
using the medium filters. We also use a pair of Chandra 
observations taken in March 2006 using the ACIS-S de- 
tector in FAINT mode, and a ROSAT observation taken 
in 1992 in pointing mode (see Table 1). 



Table 1 
Observations 



Mission 


ObsID 


Date 


Instr. 


GTI 

(8) 


GTI2 

(s) 


Count: 


exo 


7290 


2006-03-09 


ACIS-S 


27894 




300 


exo 


5436 


2006-03-11 


ACIS-S 


26799 




305 


XMM 


85280301 


2002-01-28 


MOS1 


18044 


14814 


86 








MOS2 


18051 


15070 


76 








PN 


14353 


10338 


283 


XMM 


85280801 


2002-01-30 


MOS1 


16340 


13537 


126 








MOS2 


16630 


14004 


51 








PN 


12673 


9025 


428 


ROSAT 


300181 


1992-09 


PSPCB 


45872 




452 



Note. — Observations of M13 used in this analysis, with 
GTI exposure times for XMM data quoted both with (GTI2) 
and without (GTI) aggressive background flare removal. The 
number of counts in the extraction regions are given for our 
extractions without aggressive flare removal. 

The Chandra data were reduced using Chandra In- 
teractive Analysis of Observations (CIAO) v. 4. 43 and 
Chandra Calibration Database (CALDB) v4.4.8. We 
reprocessed the data with the chandrajrepro reprocess- 
ing script to apply the latest calibration updates and 
bad pixel files, and filtered the data to the energy range 
0.3 — 10.0 keV. The Chandra data showed no strong back- 
ground flaring, so we included all data. 

The position of the qLMXB in M13 (using wavdetect 
on the combined Chandra data) is R.A.= 16 /l 41 ro 43.77 s , 
de cl.= +36°27 / 57 . 64 // , co nsistent with the position given 
in IServillat et al.1 ()2011[ ) for their Chandra source X7. 
The neighboring source, at R.A.= 16 ,t 41 m 42.47 s decl.= 
+36°28'07.29", is source X6 from IServillat et al.1 ((21)11 
(see Fig. 1), who showed it to be a CV exhibiting a 
dwarf nova eruption. Its X-ray luminosity is Lx(0.5- 
10)= 2.2^ i x 10 32 crgs/s, for a spectral shape consistent 
with an absorbed power-law of photon index F = 1.4±0.4 
(from simple spectral fits to the Chandra spectra). The 
spectra of the qLMXB were extracted using circles with 
radii of ~ 2". The specextract script generated the cor- 
responding auxiliary response files (ARFs) and redistri- 
bution matrix files (RMFs). 

The XMM-Newton data were reduced with the Sci- 
ence Analysis System (SAS) vl 1.0.0. We repipelined the 
MOS data using emchain and the PN data with epchain 
before applying the relevant filters (e.g. patterns — 12 
for MOS, - 4 for pn). We also repipelined the MOS 
an d PN data using the 2007 calibration data, as available 
to iWebb fe Barretl (|2007l ). for comparison to their anal- 
ysis (see §3.2). Both XMM observations revealed signs 
of background flaring, affecting roughly 1/3 of the obser- 
vations. Since the background flares are not extremely 
bright, we judge that spectra including all data attain a 
higher signal-to-noise ratio than if the flares were filtered 
ou t. We do, howev e r, filt er out flares to match the GTIs 
of IWebb fe Barret] ([20071 ) when reproducing their anal- 
ysis. In this case, we use filters of 5 and 10 counts per 
second for the two MOS datasets, and 30 and 50 counts 
per second for the pn data. 

Circular regions with approximate radii of 9.5" (to ex- 
clude photons from the CV X6) were used to extract 
the spectra of the NS qLMXB from the XMM observa- 

5 http:/ /cxc. cfa.harvard.edu/ciao/ 



M13 Neutron Star Atmosphere 



3 



XMM PSF HEW 


ROSAT PSF HEW 




/ '■. ". , x6 (<:v) 

/ ' 








' '/ ' \- ■■' 






M13 half-mass radius — ' 







Figure 1. Combined Chandra 0.5-6 keV image of M13, showing 
the two relevant sources X7 (the NS qLMXB) and X6 (the nearby 
CV). The radius enclosing half of M13's mass (1.69') is indicated, 
as are circles roughly indicating the half-cquivalcnt-width (HEW, 
enclosing 50% of the energy) of the point-spread functions of the 
XMM-Newton pn (radius 7.6") and ROSAT PSPC (radius 0.22') 
detectors. 

tions, except for the analysis in §3.2, where we used a 
larger (24") extraction circle, with a 13" circle around 
X6 excluded, t o matc h the extraction region used in 
IWebb fe Barretl (|2007l) . ARFs and RMFs were gener- 
ated with the SAS tasks arfgen and rmfgen. We com- 
bined the two MOS spectra in each XMM observation 
using the addspec tool to achieve better statistics, except 
for the §3.2 analysis. Finally, all spectra were grouped 
to at least 20 counts per bin. 

We used the XSELECT tool in HEASOFT0 (v6.12) 
to reduce the ROSAT data. Q. The NS spectrum was 
extracted from a circular region of approximately 47" in 
radius and then grouped to twenty counts per bin with 
no filters applied, using the on-axis response matrices. 

3. DATA ANALYSIS AND RESULTS 

We first fit the hydrogen NS atmosphere to the com- 
bined XMM, Chandra and ROSAT dataset, and calcu- 
late the constraints on the NS mass and radius. W e then 
attempt to reproduce the IWebb fe Barretl (|2007l ) result, 
using the data and calibrations available then. Finally, 
we fit helium atmosphere models, and report the NS con- 
straints in that cas e. In a l l case s, we assume a distance 
to M13 of 7.7 kpc lHarrisI () 19961 ) ( 2003 revisio n), and a 
minimum N H of 1.1 x 10 20 cm" 2 (Harris 1996). 

3.1. NSATMOS Applied to Chandra, XMM and 
ROSAT Data 

We simultaneously fit the XMM-Newton, Chandra, 
and ROSAT spectra to an a bsorbed hydrogen atmo- 
sphere NS model, NSATMOS (jHeinke et all [2006]) . The 

6 http:/ /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 

7 http: / /heasarc. gsfc.nasa.gov/docs/rosat/ros_xselect_guide/xselect_ 



NS mass was fixed to 1.4M for our initial fitting (re- 
laxed below). We use a normalization constant to allow 
for the relative difference s between each d e tector . The 
cross-calibration work of iTsujimoto et al.l ()2011l ) indi- 
cates that the MOS detectors give normalizations match- 
ing the average of several observatories, so we fix the 
MOS normalization to 1.0 and allow the pn normaliza- 
tion to float. Fixing the pn normalization to 1.0 instead 
leads to slightly (^5%) larger inferred NS radii, while fix- 
ing the Chandra normalization leads to slightly smaller 
inferred radii. 

We also included a bremsstrahlung model when fit- 
ting the ROSAT data, as ROSAT's PSPC camera was 
unable to resolve the NS a nd X6 as individual p oint 
sources (see Fig. 1), following IWebb fc Barretl (|2007l) . In 
the XMM and Chandra data sets, these sources were 
resolved, with no evidence for a contri bution by X6 to 
the N S spectrum (again agreeing with IWebb fc Barretl 
120071) . Since X6 provided few counts in the ROSAT 
data, and ROSAT has relatively poor spectral resolution 
(and no information above 2.5 keV), we fixed kT for the 
bremsstrahl ung model of X6 to 4.5 keV in our fits, con- 
sistent with IWebb fc Barretl (|2007l ) and with simple fits 
to its Chandra spectra. Omitting this component had 
relatively small effects, but our ROSAT fits include this 
component as it should be present; its fitted 0.5-10 keV 
flux was 2 x 10~ 13 ergs/cm 2 /s. 

We also tried adding a power-law component to 
the model describing all data, as found for many 
qLMXBs, with a phot on index fixed at cither 1.5 or 2 
(|Campana et al.l [19981 ). Such a power- law did not sig- 
nificantly improve the fit, and the normalization was 
consistent with zero within 90% confidence. To get the 
clearest constraints on a power- law component, we fit 
the 0.3-8 unbinned Chandra spectra with the C-statistic 
in XSPEC. This finds the power-law normalization to 
be consistent with zero, with its upper limit to be 8% 
of the total 0.5-10 keV flux, consistent with power-law 
upper limits from o ther globular cluster NS qLMXBs 
(|Heinke et al.l l2003bl ). We omit the power-law compo- 
nent from our fitting below. 

We have fit the model to both flare-filtered and un- 
filtered XMM data. While the flare-filtered data have 
less background contamination, they contain only 2/3 of 
the source photons, so the S/N ratio is similar. Table 2 
compiles results of fits to both cases, with the NS mass 
held fixed at 1.4M . We found no evidence for addi- 
tional absorption in the binary. We show a spectral fit 
to Chandra, ROSAT and unfiltcred XMM data in Fig. 2. 

We then allow the mass of the NS to vary, finding a best 
fit with NS mass 1.5M Q and radius 10.2 km (using the 
full XMM dataset). We show a contour plot of \ 2 , gener- 
ated with the steppar command, over mass and radius in 
Fig. 3. The filtered XMM dataset gave similar contours, 
though the best-fit NS mass slid down the minimal x 2 
valley to reach the lower bound of the NSATMOS model 
at 0.5M o . 

3.2. NSATMOS Applied to XMM and ROSAT Data 

We noticed that our hydrogen- atmosphere model fit, 
above, did not reac h similar results on the s ize of the M13 
NS to the work of IWebb fc Barretl (|2007l ), who require 
a relatively small NS (<11 km at 90% confidence for 
ft( m 1 afteS 1 ~1.4 M Q ). We attempted to replicate their fit by 



4 



Catuneanu et al. 



Table 2 

H Atm. Fits, CXO, XMM & ROSAT 



Parameter 



Value 



XMM Filtered for Flares 



Intrinsic Njj 
NSATMOS LogioT 
NSATMOS R 
Reduced x 2 /dof 
Null Hyp. Prob. 



0.0 



+0 l X 10 20 cm 
+0.05 



5-99%.03 
+ 19 km 



11.7 



0.7551/73 
0.94 



Without Flare Filtering 



Intrinsic Nh 
NSATMOS Log 10 T 
NSATMOS R 
Reduced \- 2 /dof 
Null Prob. 



0.0 



+0.7 

0.0 

6.00 



x 10 2 



+0.06 
0.04 



10.6 



+2.1 



km 



0.8337/84 
0.86 



Note. - Chandra, ROSAT and XMM-Newton (flare- 
filtered and unfiltered) data fit to an NSATMOS hydrogen- 
atmosphere NS model. A bremsstrahlung model was added 
to the ROSAT spectrum to model X6's spectrum, included in 
the ROSAT extraction. The second fit corresponds to Fig. 3 
for a NS mass fixed to 1.4M0. The Nh quoted is any intrinsic 
Nh hi the binary, in addition to the (fixed) cluster value of 
1.1 X 10 20 cm" 2 . 
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Figure 3. Mass vs radius confidence contours for the non-flarc- 
filtcrcd fit to XMM, Chandra and ROSAT data with the NSAT- 
MOS hydrogen atmosphere. The solid (blue), dotted (green), 
and dashed (red) lines represent 90%, 99%, and 68% confidence 
contours respec tively. The upper left section (with R < 2.8M, 
Lindblom 1984) is shaded to indicate its inaccessibility for any NS, 
as any equation of state would require the sound speed to exceed 
c. 
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Figure 2. Spectra of M13 qLMXB, including unfiltered XMM 
MOS (red and black) and pn (blue and green), Chandra (light 
blue and magenta), and ROSAT (light gray) data, fit with the 
NSATMOS hydrogen-atmosphere model. (The ROSAT model also 
contains a bremsstrahlung component, for the unresolved CV X6.) 
The lower portion of the diagram depicts x 2 residuals. 

reprocessing and extracting the XMM data as in Webb 
et al. (see §2 above). For a NS mass fixed to 1.4 M Q , we 
obtain the fit in Table 3. 

When allowing the mass to vary, we find a best-fit NS 
mass of 1.3M© and a 10.0 km radius, not dissimilar to 
the best fit of Webb et al. (their Table 2). However, 
our contours (Fig. 4) are rather less constraining than 
those in Webb et al.'s Fig. 6. Webb et al.'s quoted best 
fits lie very close to their upper 90% confidence contour 
lines. We do not see such behavior in any well-behaved 
chi-squared contour plots of the M13 qLMXB, or any 
other NSs in our experience. The uncertainties on radius 
quoted in their Table 2 also seem unusually narrow. 



Table 3 

H Atm. Fits, XMM & ROSAT 



Parameter 



Value 



Intrinsic Nh 
NSATMOS LogioT 
NSATMOS R 
Reduced x 2 /dcf 
Null Hyp. Prob. 



0.0 



-0.6 
-0.0 



x 10 



20 



i.o r 



+0.06 
0.07 

+ i 5. 
0.8595/68 
0.79 



9.5 



Note. — ROSAT and flare-filtered XMM spectra 
of the M13 NS fit by an NSATMOS model, with a 
bremsstrahlung component included to fit X6 in the 
ROSAT data. This corresponds to Fig. 4 for a NS 
mass fixed to 1.4M0. Errors marked with * indicate 
the parameter hits a hard limit of the model. The Nh 
quoted is any intrinsic Nh in the binary, in addition 
to the (fixed) cluster value of 1.1 X 10 20 cm" 2 . 



3.3. A Helium Atmosphere 

Finally, we consider the effects of a helium atmosphere 
on the inferred NS propert ies. We used the hel ium atmo- 
sphere model described in IHo fc Heinkel (|2009|) to fit the 
data as described in §3.1. We present fits with the NS 
mass fixed to 1.4 M , using flare-filtered vs. unfiltered 
XMM data, in Table 4, and plot mass vs. radius contours 
in Fig. 5. 

Although the flare-filtered data have a slightly better 
fit to the hydrogen model {8\ 2 = 1.5 for the same degrees 
of freedom) , the helium model fit is still acceptable. The 
fits to the unfiltered XMM data produce essentially indis- 
tinguishable x 2 values for helium or hydrogen fits. Thus, 
we cannot determine from the current data whether a 
hydrogen or helium atmosphere is the correct one, and 
thus which mass-radius contours arc appropriate for the 
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Figure 4. Mass vs. radius contours from the NSATMOS 
hydrog en-atmosphere fi t to XM M and ROSAT data processed fol- 
lowing IWebtTEBarret] 420071) . The solid (blue), dotted (green) 
and dashed (red) lines represent 90%, 99%, and 68% confidence 
contours respectively. 



Table 4 

Helium Atm. Fits, CXO, XMM & ROSAT 



Parameter 



Value 



XMM Filtered for Flares 



Intrinsic Njj 
He NS LogioT 
He NS Ft 
Reduced x 2 /dof 
Null Prob. 



5.94 



x 10 20 cm" 

0.05 
0.04 



14.6 



+3.5 



km 



3.1 

0.7753/73 
0.92 



Without Flare Filtering 



Intrinsic Nh 
He NS LogioT 
He NS R 
Reduced x 2 /dof 
Null Prob. 



0.2 



+0 i X 10 20 en 

5 96+ 05 
°- yD -0.05 

I2.8ti; 2 km 
0.8341/84 
0.86 



Note. — Fits to our helium atmosphere model. The second 
fit corresponds to Fig. 5 for a NS mass fixed to 1.4 MQ(see 
text). The Nh quoted is the intrinsic Nh in the binary, added 
to the (fixed) cluster value of 1.1 X 10 20 cm -2 . 



M13 qLMXB. 



4. DISCUSSION 



Whether a NS qLMXB hosts a hydrogen or helium 
atmosphere affects the mass and radius constraints in- 
ferred from fitting its quiescent X-ray spectrum. For 
the M13 qLMXB, for instance, the best-fit NS radius in- 
creases by ^2 km with a helium atmosphere. Although 
ultracompact LMXBs arc known to be common in glob- 
ular clusters, it is not certain whether ultracompact 
qLMXBs will display pure helium (or, for hybrid w hite 
dwarf donors, carbon; see lNelemans fc Jo nker 2010fl at- 
mospheres, which introduces uncertainty into calcula- 
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Figure 5. Mass vs. radius contours when fitting the Chandra, 
non-flarc-filtcrcd XMM, and ROSAT data with a helium atmo- 
sphere model. 

tions of the radii of any NSs that may be in ultracompact 
qLMXBs. 

There are several ways to distinguish ultracompact 
qLMXBs from normal qLMXBs. Detection of hydrogen 
in the spectrum, or (narrow-band) photometry, of the 
companion can provide a clear distinction; this proves 
that the qLMXB in uj Ce n has a hydrogen atmosphere 
(e.g. lHaggard et al.|[2004T ). Measurement of the orbital 
period (e.g. through eclipses or pulsation timing in 
outburst) allows clear discriminatio n, as done for the 
eclipse rs W37 and X5 in 47 Tuc (jHeinke et al.l l200l 
I2003al) and the millisecond pul sars SAX J1748.9-2021 
and N GC 6440 X-2 in NGC 6440 (jAltamirano et al.|[2008l 
l2010f) . The energetics and duration of X-ray bursts, when 
these can be clearly attributed to a particular qLMXB, 
may provide evidence for or against the presence of hy- 
drogen; e.g., the surface of the Terzan 5 transient EXO 
1745- 248 is known to contain hydrogen (jGallowav et al.l 
120081 ). For the M13 qLMXB, narrow-band photometry 
should be the first project. Finding evidence of H-a emis- 
sion from an optical counterpart, using the Chandra po- 
sition and archival (or new, deeper) Hubble Space Tele- 
scope observations of M13, would prove that the NS pho- 
tosphere is made of hydrogen. 

Another concern is the presence of cross-calibration 
differences between the various detectors used for this 
type of work; Chand ra's ACIS detector, and XMM's pn 
and MOS detectors. iTsuiimoto et al.l (|2Qlll ) found that 
the XMM pn detector gives 1-8 keV fluxes 6.5% lower 
than the MOS detectors in simultaneous observations, 
while Chandra's ACIS-S3 detector averages 11.6% higher 
fluxes. Our analysis gave average XMM pn normaliza- 
tions 5.5% higher, and average Chandra ACIS-S3 nor- 
malizations 12.5% higher, than the MOS normalizations, 
consistent with other Chandra /XMM cross-c alibration in 
the 0.5-2 keV range (jNevalainen et al.l[2010| ). Fixing the 
other instruments, instead, to have normalization equal 
to one would systematically increase the inferred NS ra- 
dius by roughly 3 and 6%, respectively. It is not obvious 
which detector's calibration is more accurate, which sug- 
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gests a ^5% systematic uncertainty in all such NS radius 
measurements . 

Our resul t affects inferences of th e equation of state of 
NSs, as the iWebb fc Barret! (|2007t ) analysis of the M13 
NS gave some of the tightest constraints, consistent with 
the IQzel et all (|2010h meta-analysis of X-ray burst re- 
sults that preferred a relatively small NS radius under 
10 km. Our re-analysis does not support small radii for 
a hydrogen atmosphere fit. For a helium atmosphere, 
the radius would be even larger, >10 km at 1 sigma for 
typical NS masses. Our results are consistent with the 
ISteiner et al.l (|2010l ) meta-analysis of NS bursters and 
qLMXBs, which pr eferred a NS radius between 11-12 km. 
ISteiner et al.l 120121 explicitly sho w that their results are 
robust against the removal of the IWebb fc Barretl (|2007l ) 
M13 qLMXB constraints. Additional high-quality NS 
constraints would be very valuable in further constrain- 
ing the NS equation of state, but should consider un- 
certainties in atmospheric composition and absolute flux 
calibration. 
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